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Abstract 

The cross section of the diffractive Drell-Yan (DY) process, pp — > 
UXp, where the system UX is separated by a large rapidity gap from 
the recoil proton, is calculated in the light-cone dipole approach. This 
process reveals unusual features, quite different from what is known 
for diffractive DIS and nonabelian radiation: (i) the diffractive radi- 
ation of a heavy dilepton by a quark vanishes in the forward direc- 
tion; (ii) the diffractive production of a dilepton is controlled by the 
large hadronic radius; (iii) in contrast with DIS where diffraction is 
predominantly soft, the diffractive DY reaction is semihard-semisoft; 
(iv) as a result of the saturated shape of the dipole cross section, the 
fraction of diffractive DY events steeply falls with energy but rises 
as function of the hard scale. These features are common for other 
abelian bremsstrahlung processes (higgsstrahlung, Z-strahlung, etc.). 
Measurements of diffractive DY processes at modern colliders would 
be a sensitive probe for the shape of the dipole cross section at large 
separations. 

1 Diffractive radiation: heuristic approach 

Diffraction excitation of hadrons is possible due to the presence of quan- 
tum fluctuations in the projectile particles. In classical physics only elastic 



1 



diffraction is possible. It was first realized by Feinberg and Pomeranchuk pQ 
and Good and Walker [2J that the compositeness of hadrons leads to pro- 
duction of new states. Although different Fock components of the hadron 
experience only elastic scattering, which is a shadow of inelastic collisions, 
the wave packet composition may be altered producing a new hadronic state. 
Indeed, this may happen if the Fock states interact differently, otherwise the 
wave packet retains the same composition, i.e. the final and initial states are 
identical. 

The dipole description of diffraction in QCD was presented in [3l S] - Since 
dipoles of different transverse size tt interact with different cross sections 
cr(r<r), this gives rise to single inelastic diffraction with a cross section given 
by the dispersion of the r T -distribution [3J, 



Here is the transverse momentum of the recoil proton; (<7(ty)) is the 
dipole-proton cross section averaged over dipole separation. 

The dipole description of diffractive radiation of photons and gluons was 
developed in Ref. [5]. Diffractive radiation of a photon by a quark, in which 
the photon can be either real, or heavy decaying into a dilepton, turns out to 
vanish in the forward direction. Indeed, let us consider two Fock components 
of a quark, just a bare quark, \q), and a quark accompanied by a Weizsacker- 
Williams photon, [qj*). In both components only the quark can interact, 
therefore the two terms in ([1]) cancel each other. Notice that the partial 
diffractive amplitude at a given impact parameter does not vanish, since the 
recoil quark in the \q^y*) state gets a shift in impact parameters compared to 
the \q) state, and the two Fock components interact differently. Only after 
integration over impact parameter, corresponding to the forward amplitude, 
the diffractive radiation vanishes. 

A direct calculation of Feynman graphs [5] confirms this expectation. 
Notice that even nondiffr active photon radiation in inelastic collisions is im- 
possible without momentum transfer. Indeed, the Born graphs Fig. dh and 
[Tb cancel if px = 0. Intuitively it is clear that if the electric charge gets no 
kick and is not accelerated, then no radiation happens. 

This result is retained in the nonabelian case, namely a quark does not 
radiate gluons if the t-channel gluon provides no momentum transfer [6]. 
This is not a trivial result, since in this case a color current flows between 
the beam and target. 
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Figure 1: Born approximation for radiation of a photon by a quark in inelastic 
collision 



In the case of diffraction forward electromagnetic radiation vanishes as 
well, but this is less obvious. Indeed, even if the two-gluon exchange provides 
no momentum transfer, each of gluons can carry a transverse momentum. 
The relevant Born graphs for diffractive photon radiation are shown in Fig. [2j 
Graph (b) does not contribute to radiation, provided that the radiation time 
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Figure 2: Born graphs for diffractive electromagnetic radiation by a quark. 



considerably exceeds the duration time of interaction. Only graphs (a) and 
(c) can contribute, but they cause no radiation for forward scattering for the 
same reason as in the inelastic collision, Fig. [TJ 

This can be interpreted as a consequence of Landau-Pomeranchuk prin- 
ciple [7j : radiation depends on the whole strength of the kick, rather than on 
its structure, if the time scale of the kick is shorter than the radiation time. 
In other words, if two opposite kicks are separated by a short time interval, 
the radiation spectra from each of the two kicks interfere destructively and 
compensate each other [8j[9], i.e. no radiation occurs. 

Notice that disappearance of abelian diffractive radiation in the forward 
direction goes along with the result of [10J that the diffractive cross section is 
proportional to the mean momentum transfer squared. That was, however, 
an oversimplified treatment of the Pomeron as a point-like vector meson. 

Although a quark cannot radiate a dilepton in forward direction, a hadron 
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can. That is possible due to transverse motion of the valence quarks in the 
hadron, i.e. abelian radiation even at a hard scale is sensitive to the hadron 
size, which is a dramatic break down of QCD factorization [11] (which has 
never been proven for this process). Failure of factorization for diffractive 
Drell-Yan reaction has been already known. It was found in [T2], [13] that 
factorization fails due to the presence in the Pomeron of spectator partons. 
Below we demonstrate that factorization in Drell-Yan diffraction is even more 
broken due to presence of spectator partons in the colliding hadrons. 

Notice that diffractive nonabelian (gluon) radiation is different, it does 
not vanish in the forward direction [5J. This is a direct manifestation of the 
nonabelian dynamics. In addition to the three graphs for electromagnetic 
radiation shown in Fig. [21 in the case of gluon radiation there is a fourth 
term corresponding to one of the t-channel gluons coupled to the radiated 
one. This term gives rise to a nonzero forward diffraction. This is confirmed 
by data, since a nonzero forward cross section of diffractive gluon radiation, 
called triple- Pomeron term in Regge approach, is well established experimen- 
tally [33]. Nevertheless, the cross section of diffractive gluon radiation turns 
out to be amazingly small. Indeed the Pomeron-proton total cross section 
extracted from data for large mass diffraction turns out to be only 2 mb, an 
order of magnitude smaller than for pion-proton. This effect is discussed and 
explained in Ref. [5]. 



2 Calculation of Feynman graphs 

Although a quark cannot diffractively radiate photons in forward direction, 
a hadron can. Let us consider a pp collision, 

Pb + Pt -> 7* + X + p t , (2) 

assuming that the initial beam nucleon and the final state consist of three 
valence quarks, pb = 3^; X = 3g^. The amplitude of the Drell-Yan (DY) 
reaction Eq. (j2j) has three terms, 

Af = 4§ + 4? + Aff , (3) 

where subscripts i and / correspond to initial and final states of the Drell-Yan 
diffractive reaction, Eq. (J2J). 

Each term in ([3]) corresponds to radiation of the 7* by one of the valence 
quarks, and it includes contributions from tree graphs shown for A^) in 
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Fig. [31 The first graph (a) leads to no radiation at px — as was explained 
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Figure 3: Born graphs for diffractive electromagnetic radiation by a proton 
above. The rest, graphs (b) and (c), can be calculated according to, 
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Here are the light-cone wave functions of the 3q systems in the initial and 
final state respectively; r 1; r* 2 and r 3 are the impact parameters of the quarks; 
x qi ,x q2 ,x q3 are the fractions of the proton light-cone momentum carried by 
the quarks; k and the photon transverse and fractional longitudinal 

momenta; r is the transverse separation between the photon and the radiating 
quark; a = x 1 jx qx \ and $j(r, a) is the distribution amplitude for photon 
radiation by quark g« in the mixed representation, transverse coordinate and 
longitudinal momentum. The factor has the form, 



ri,r 2 ,r 3 ,r,a) 



a{r 1 -r 2 



■a{f\ — r 2 — ar)+a(fi — 7%) — cr(fi — r 3 — ar) , 

(5) 

where airr) is the universal dipole cross section [3]. 

The structure of this amplitude is easy to understand. According to 
the discussion in the previous section, the amplitude of diffractive radiation 
should be proportional to the difference between elastic amplitudes of the 
two Fock components, with and without the photon. In both cases only the 
3-quark dipoles interact, but they have different sizes in the |3g) and |3g7*) 
components, therefore interact differently. Indeed, when a quark fluctuates 
into quark-photon with transverse separation r, the final quark gets a trans- 
verse shift Ar = ar, where a is the fraction of the quark momentum taken 
away by the photon. This shift is explicitly presented in ([S]). 
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For the model of symmetric proton wave function we use, the interference 
terms in the cross section cancel, since they are proportional to Z\Z 2 + ZiZ^ + 
Z2Z3 = (no cancellation occurs for neutrons), where Zj is the electric charge 
of the ith quark. Thus, summing over different final states / one gets for the 
forward diffractive cross section, 
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Employing completeness, 
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we arrive at 
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Here $(r, a) = <&j(f, a)/Zj is the amplitude of photon radiation by a quark 
with charge one. 

To progress further we assume the Gaussian shape for the quark distri- 
bution in the proton, 
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x (J(fi + f 2 + f 3 ) 5(1 - - x 92 - i B ) , (9) 



where a = (r^) -1 is the inverse proton mean charge radius squared. The 
three-body distribution function p(x qi , x q2 ,x q3 ) should reproduce the valence 
quark distribution in the proton, 



J dx q2 dXq 3 p(x qi , Xg 2 , X q3 ) p qi (x qi ) 



(10) 



Summing different quark and antiquark species one arrives at the proton 
structure function [15J, 

£z£ [ Pq ( x )+ Pq ( x )} = l-F 2 (x) . (11) 

q X 

At this point we generalize our three-body quark wave function to a 
multi-body one including antiquarks. However, we keep the same simple 
coordinate part of the wave function. 

In what follows we rely on the popular saturated form of the dipole cross 
section, 

a{r) = a (l - e 1 *'**) (12) 

where the parameters fitted to DIS data at small x can be found in [16J. 
The distribution functions have the form, 
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for radiation of longitudinally and transversely polarized photons respec- 
tively. Here K (rr) is the modified Bessel function; 
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Xij are the spinors assigned to the quark before and after the radiation 
respectively; a is the Pauli matrix; n is a unit vector directed along the 
initial collision momentum; m q is the quark mass. 

Now we are in a position to perform integration over coordinates in (jSJ). 
Using the integral representation for the modified Bessel function, 
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we arrive at the longitudinal cross section of diffractive DY reaction, 
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The cross section of radiation of transversely polarized dileptons summed 
over polarizations reads 
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Here Xi = x 7 and x 2 are the standard DY variables, 
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s s 
where q and P\^ 2 are the 4-momenta of the radiated dilepton, and projec- 
tile/target protons respectively. Other notations in (ITS]) are, 
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and i = 1,2, 3, 4. 

Since the forward cross section is known, the total diffractive cross section 
can be estimated as, 
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DY 
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(21) 
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where B D /(s) is the slope of the t- dependence of the cross section, which is 
similar to the corresponding slope measured in diffractive DIS. 

3 Unitarity corrections 

While the radiated leptons do not interact with the target, the accompanying 
quarks do, and they can easily initiate an inelastic interaction followed by 
multiparticle production filling up the rapidity gap. Therefore, the survival 
probability of the gap causes a suppression for the diffractive cross section. 
This is easy to take into account in the impact parameter representation. 
The amplitude Eq. ([3]) acquires a factor, 

A lf (b) ^ A lf (b) [1-ImOT , (22) 

where f^(b) is the partial elastic amplitude. After squaring this amplitude 
and integration over impact parameter we arrive at the cross section which 
is different from one in Eq. (12~T1) by a suppression factor K, do®/ ' jdx\ =>- 
K do® J I dxi, where [17] 

1 n B®/(s) + 2Blf{s) + (4^ B™{s) [B®/(s) + B^(s)} 

Here the elastic slope depends on energy as B^f(s) = B® t + 2a^p ln(s/s ) 
with B% = 7.5 GeV" 2 , s = 1 GeV 2 . The slope of single-diffractive DY cross 
section can be estimated as, B D /(s) ~ ( r ch)/3 + 2ajp ln(s/so), where the 
proton mean charge radius squared (r 2 h ) = 0.8 fm 2 . 

Notice that this eikonal expression for the gap survival probability is a 
conservative estimate. Inclusion of corrections related to intermediate diffrac- 
tive excitations of the beam particle (so called Gribov corrections [18]) make 
the medium more transparent, i.e. increase the survival probability of the 
rapidity gap [T7] . 

Our estimate Eq. fl2U|) can be compared with the results of more elaborate 
models [T^ [201 ED E2] incorporating a part of the Gribov corrections. The 
predicted suppression factors have similar orders of magnitude. Notice that 
one can easily replace the suppression factor (|23|) by a preferable one. 
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4 Diffractive vs Inclusive 



We should compare the diffractive DY cross section with the inclusive one, 
which also can be calculated within the dipole approach [23j I2U EES] , 
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This simple formula reproduces with high precision the cross section mea- 
sured by E772 [25] and E866 [26J experiments at Fermilab, and agrees quite 
well with the results obtained from much more complicated NLO parton 
model calculations [27 . 

Applying Eqs. to Eq. (EP we get, 
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The notations here are similar to ([2"01 . 

o = (l- y 2 )(2w + l) ; 
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We perform our calculations with the proton structure function parametrized 



as, 



with Ql 



F 2 (x,Q 2 ) = A(x) 
20 GeV 2 and A = 



ln(Q 2 /A 2 ) 
ln(Q 2 /A 2 ) 



fl(x) 



1 + 



C{x) 

Q 2 



(28) 



0.25 GeV. This parametrization was fitted 
to available data for F 2 (x,Q 2 ) from different experiments. The form and 
fitted parameters for the functions A(x), B(x) and C(x) can be found in the 
Appendix of [28] . 
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For the parameters of the dipole cross section, Eq. (fT21) . we use the results 
of the fit of Ref. [B] to HERA data: a = 23.03 fm; R (x 2 ) = 0.4 fm x 
(x 2 /x )°- U4 , where x = 3.04 x 10" 4 . 

We calculated the ratio of the diffractive to inclusive cross sections. The 
results are plotted in Fig. 0] as function of the dilepton effective mass squared 
at fixed x l = 0.5, 0.9 and energies y/s = 40 GeV, 500 GeV and 14TeV. In 
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Figure 4: Ratio of diffractive to inclusive DY cross sections as function of 
M 2 at fixed x\ = 0.5 (solid curves) and x\ = 0.9 (dashed) and energies (from 
top to bottom) y/s = 40 GeV, 500 GeV and 14 TeV. 

these calculations we summed the contributions of longitudinal and trans- 
verse parts both in the diffractive and inclusive cross sections. 

The results depicted in Fig. H] demonstrate unusual effects. The diffractive- 
to-inclusive cross section ratio is steeply falling with energy. This is counter- 
intuitive, since diffraction, which is proportional to the dipole cross section 
squared, should rise with energy steeper than the total inclusive cross section. 
At the same time, the ratio rises with the hard scale, M 2 . This also looks 
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strange, since diffraction is usually associated with soft interactions 

To understand these remarkable features of DY diffraction we should 
return back to the original diffractive amplitude, Eq. (j!])-©. As we empha- 
sized, the diffractive amplitude is given by the difference between the cross 
sections of the relevant Fock states with and without the dilepton. This is 
explicitly incorporated into Eq. (jSJ). Assuming r ~ 1/M <C Rq(x2) we can 
expand this cross section difference as, 

a(R) - a(R - af) = e -« 2 /^ 2 ) ff . R) + 0(r 2 ) . (29) 

This is an interesting result: the amplitude is linear in r, i.e. the diffractive 
cross section is quadratic function of r. This is different from diffraction in 
DIS where the cross section is oc r 4 . The latter is predominantly a soft pro- 
cess, since the end-point qq fluctuations, a — > 0, 1, have no scale dependence, 
only their weight is oc 1/Q 2 - Therefore these soft fluctuations dominate the 
diffractive DIS cross section [291 EQ]- However, since diffractive DY cross 
section is oc r 2 , soft and hard interactions contribute on the same footing, 
and their interplay does not depend on the scale, similar to inclusive DY or 
DIS. 

Moreover, expression fl30l) leads in Eq. (jSJ) to the same integral over r and 
a as in the inclusive cross section Eq. (|24|) . Therefore, all energy and scale 
dependence of the diffractive-to-inclusive cross section ratio comes via the 
x 2 -dependent factor, 

rr DY 1 

°% RIM ' ( j 

Since for light hadrons ^0(^2) < 2-R 2 , this expression rises with Rq, i.e. with 
X2 = M 2 jx\S. This explains why the ratio depicted in Fig. H] falls down 
with energy, but rises with M 2 . Note that the falling energy dependence is 
partially due to the absorptive corrections (j23|) . 

We found that both the diffractive and inclusive DY cross sections are 
dominated by radiation of transversely polarized DY pairs. This is demon- 
strated in Fig. [5] where we plotted the ratios of the longitudinal to transverse 
cross sections for inclusive diffractive DY reaction as function of M 2 at dif- 
ferent energies. These ratios turn out to be very similar to that for inclusive 
DY processes. 

We see that the diffractive ratios are close to the inclusive ones. They do 
not exhibit any strong dependence on energy, but fall down at harder scales, 
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Figure 5: Ratio of longitudinal to transverse parts of diffractive DY cross 
sections as function of M 2 at fixed x\ = 0.5. The predictions are made for 
energies, y/s = 40 GeV (dotted), y/s = 500 GeV (dashed), and y/a = 14TeV 
(solid). 

similar to DIS [31J. This fall is related to the suppression of the end-point, 
a — > 1, of longitudinally polarized DY fluctuations in equation (113 1) . 

5 Summary 

Contrary to the simple intuition based on QCD factorization, diffractive DY 
processes have properties quite different from what is known about DIS and 
nonabelian radiation. 

• A quark cannot radiate noninteracting particles (photons, dileptons, 
gauge bosons, Higgs, etc.) diffractively in the forward direction, i.e. 
without gaining any momentum transfer. This is a manifestation of the 
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general principle: if a charge is accelerated and then immediately decel- 
erated, these two sources of radiation interfere destructively and cancel 
each other [81 |9] . This can also be interpreted in terms of the Landau- 
Pomeranchuk principle: radiation at times longer than the time interval 
of the interaction depends only on the strength of the whole kick, but 
does not resolve its structure (a single or multiple kicks). 

The nonabelian case, QCD, is different: a quark can radiate gluons 
diffractively in the forward direction^- This happens due to possibility 
of interaction between the radiated gluon and the target. 

• Nevertheless, a hadron can radiate abelian fields diffractively and in 
the forward direction. This process involves the spectator quarks as 
is illustrated in Fig. [3X and for this reason the cross section depends 
on the hadronic size and rises with it. This is an apparent and strong 
breakdown of factorization. 

The situation in the nonabelian case is different. The main contribution 
to the cross section comes from diffractive gluon radiation by a single 
quark. Interaction with the spectator quarks has little impact and 
vanishes at large separations. 

• The ratio of the cross section of diffractive radiation of abelian fields 
to the inclusive one is steeply falling with energy as depicted in Fig. HI 
This occurs even without unitarity corrections which make the fall even 
steeper, and is a result of the saturated shape of the dipole cross section. 
If this cross section was rising like <r(r) oc r 2 , the ratio would be energy 
independent (up to unitarity corrections). This property is strikingly 
opposite to what would follow from factorization [TT] . 

This again is different in the nonabelian case, since here without unitar- 
ity corrections the diffractive cross section rises faster than the inclusive 
one. This happens because the diffraction cross section is a quadratic 
function of the dipole cross section, while the inclusive is linear . 

This is also different from DIS where the diffraction to inclusive ratio 
rises with energy. 

1 This is a higher order effect. In Born approximation for inelastic collisions the forward 
gluon radiation vanishes [6] . 
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• Unexpectedly, the fraction of diffractive events to the total inclusive 
cross section of abelian radiation rises with the dilepton effective mass 
M (see Fig. H]). This effect has the same origin as the energy depen- 
dence, the dipole cross section which levels off at large separations. On 
the contrary, the fraction of diffraction in the total DIS cross section 
is slowly, logarithmically, falling with the hard scale. Same would be 
true for Drell-Yan diffraction if factorization were correct. 

The rise of diffraction with the hard scale looks counterintuitive, since 
diffraction is usually associated with soft interactions. However, that 
is not true for diffractive DY processes. 

• Hard and soft interactions contribute to DY diffraction on the same 
footing, and their ratio is scale independent like in inclusive DY process. 
This is a result of the specific property of DY diffraction: its cross 
section is a linear, rather than quadratic function of the dipole cross 
section. 

On the contrary, diffractive DIS is predominantly a soft process, be- 
cause its cross section is proportional to a 2 (r) 

• The main features of our results for Drell-Yan reaction are valid for 
other diffractive abelian processes, like production of direct photons, 
Higgsstrahlung, radiation of Z and W bosons. One should not rely on 
QCD factorization for these reactions. 

• Measurement of diffractive DY processes in a wide energy range from 
fixed target experiments up to the modern colliders, RHIC, Tevatron 
and LHC, would provide a precious direct probe for the behavior of 
the dipole cross section at large distances, above the saturated radius, 
r > Rq. 
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